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The mounting evidence for anomalously large peculiar velocities in our Universe presents a chal- 
lenge for the ACDM paradigm. The recent estimates of the large scale bulk flow by Watkins et al. 
are inconsistent at the nearly 3cr level with ACDM predictions. Meanwhile, Lee and Komatsu have 
recently estimated that the occurrence of high- velocity merging systems such as the Bullet Cluster 
(1E0657-57) is unlikely at a 6.5 — 5.8<r level, with an estimated probability between 3.3 x 10 -11 and 
3.6 x 10 -9 in ACDM cosmology. We show that these anomalies are alleviated in a broad class of 
infrared-modifed gravity theories, called brane-induced gravity, in which gravity becomes higher- 
dimensional at ultra large distances. These theories include additional scalar forces that enhance 
gravitational attraction and therefore speed up structure formation at late times and on sufficiently 
large scales. The peculiar velocities are enhanced by 24—34% compared to standard gravity, with the 
maximal enhancement nearly consistent at the 2a level with bulk flow observations. The occurrence 
of the Bullet Cluster in these theories is « 10 4 times more probable than in ACDM cosmology. 



A striking discrepancy with ACDM predictions has 
emerged in recent estimates of the local large scale bulk 
flow. Using a compilation of peculiar velocity surveys, 
Watkins et al. [1 , 2 finds a bulk flow of 407 ± 81 km s" 1 
on 50 h~ 1 Mpc scales, which is inconsistent at the ps 3<t 
level with the ACDM rms expectation of ^180 km/s. 

Evidence for this anomaly on larger scales, though 
less reliable, comes from measurements of cluster pecu- 
liar velocities using the kinetic Sunyaev-Zel'dovich ef- 
fect. These indicate a coherent bulk motion of 600- 
1000 km s' 1 out to > 300/i -1 Mpc 011]. These results 
match the direction of the flow found by pQ on the length 
scales where they overlap. See also [5] for recent efforts 
in reconstructing the local group peculiar velocity. In 
the future, observations beyond our local area will offer 
better statistics on bulk flows [5]. 

Independent evidence that structure is evolving more 
rapidly than expected is the "Bullet Cluster" 1E0657- 
57 [7J. A key input in hydrodynamical simulations of 
this system [8HTU] is the initial subcluster velocity. Re- 
cent simulations have shown that an initial velocity of 
3000 km/s at 5 Mpc separation is required to explain 
the data (TT]. Using horizon scale N-body simulations, 
Lee and Komatsu |T2] have estimated that the probabil- 
ity of having such large velocities in ACDM cosmology 
is between 3.3 x 10~ n and 3.6 x 10 -9 — that is, the 
"Bullet" system is between 6.5 and 5.8<r away from the 
mean velocity for colliding clusters. A previous estimate 
of this likelihood has been interpreted as evidence for a 
new attractive force in the dark sector [13] . See [21 [15] 
for examples of other violent merging systems. 

These anomalies motivate us to study peculiar veloc- 
ities in a broad class of infrared (IR) modified gravity 
theories called brane-induced gravity [L6H26] . See [27] for 
a useful parametrization of modified peculiar velocities, 
and [55] for related work. The theories of interest involve 
extra scalar degrees of freedom in AD, which are inherited 
from the higher-dimensional massless graviton. These de- 
grees of freedom couple to the trace of the matter stress- 



energy tensor and therefore enhance the effective gravita- 
tional attraction compared to Newtonian gravity at late 
times and on large scales. Large scale structure is more 
developed and is presently evolving faster on large scales 
than in the ACDM model, both of which lead to larger 
bulk flows [291 US- 

While we focus on higher-dimensional theories for con- 
creteness, enhanced peculiar velocities are also expected 
in any theory with a long range extra scalar force, such 
as galileon scalar-tensor theories [31TI37] , interacting dark 
sector models [55H14] and symmetron theories [35J 05] . 
In chameleon// (R) [4TH55] models, however, the scalar 
fifth force has a range $ Mpc and hence cannot explain 
the large scale anomalies discussed here. It would be in- 
teresting to study whether our results also apply to the 
recent IR- modified theory proposed by [5"4"] . 



A. Brane-induced Gravity 

The principal motivation for modifying General Rel- 
ativity (GR) at ultra-large distances is the cosmologi- 
cal constant problem |55j . (See [56) for a recent review 
of cosmological tests of gravity.) Vacuum energy is the 
zero-momentum component of stress energy and hence its 
backreaction depends sensitively on the nature of gravity 
in the far infrared. A compelling approach is degravita- 
tion, in which gravity acts as a high-pass filter [5TH59"] : 
the cosmological term is in fact large, in accordance with 
field theory expectations, but gravitates very weakly. 

A perennial challenge in devising consistent IR- 
modified theories of gravity is quantum stability, i.e. 
avoidance of ghost-like (negative energy) instabilities. 
Giving the graviton a hard mass a la Pauli-Fiertz [60] . 
for instance, unavoidably leads to instabilities [6TM63| . A 
more promising approach is brane-induced gravity [161 - 
l26l I64H661 l69l [70] , which relies on branes and extra di- 
mensions. The most widely known example is the Dvali- 
Gabadadze-Porrati (DGP) brane- world model [ToTll8| . 
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The normal branch of the DGP model is perturbatively 
ghost free, in contrast to the self-accelerating branch [7H - 
[75], and thus represents a perturbatively consistent IR 
modification of gravity. 

The Cascading Gravity framework [2"TH2"4l R)9"] extends 
the DGP model to D > 6 space-time dimensions. In 
the simplest version with a 6D bulk space-time, our 3- 
brane is embedded in a 4-brane, each with their own 
induced gravity terms. The upshot of this generaliza- 
tion is twofold. First, the soft mass term for the gravi- 
ton is a more slowly varying function of momentum than 
in DGP, which is a necessary condition for degravita- 
tion |59| . Thus Cascading Gravity is a promising frame- 
work for realizing this phenomenon [23] . Furthermore, 
the cascading graviton mass term results in an expansion 
history that closely resembles ACDM cosmology and is 
therefore less constrained by observations [291 El] ■ 

What about perturbative stability? Perturbing around 
6D Minkowski space with empty branes reveals a ghost 
scalar mode. Early work [21] revealed, however, that 
the ghost is excised by including a sufficiently large 
tension on the 3-brane or, alternatively, by consider- 
ing a higher-dimensional Einstcin-Hilbert term localized 
on the brane [22, 66 -68 . While the original derivation 
of [2T] was restricted to a particular decoupling limit of 
the theory, recently the absence of perturbative ghosts 
has been proven rigorously by perturbing the full 6D so- 
lution in the presence of brane tension 24J . These results 
establish the cascading gravity framework as a perturba- 
tively consistent IR modification of gravity. 



B. Summary of Results 

In this paper we show that the observational bulk flow 
anomaly is alleviated in brane-induced gravity theories. 
The bulk flow enhancement depends on the number D 
of bulk space-time dimensions and the cross-over scale 
r c beyond which gravity on the brane becomes higher- 
dimensional. Our bulk flow results can be summarized in 
a fitting formula for the rms of the 1-dimensional velocity 
on 50ft -1 Mpc scales, which is valid for D > 4 and r c < 
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The superscript "G" indicates the use of a Gaussian win- 
dow function. In particular, for D = 4 (i.e., 7 = 2/3), 
this matches the ACDM prediction: = 102 km/s 

for our fiducial parameter choices. In Sec. Ill B| we gen- 
eralize ( B to a fitting formula valid on a range of scales 
- see (23). Our fiducial cosmology is consistent with 



WMAP 7-year data [75] and assumes a spatially-flat uni- 
verse with tt m = 0.24, Ob = 0.042, a primordial power 
spectrum with tilt n s — 0.96, and a primordial amplitude 
chosen to yield a present-day normalization of ag = 0.8 
for a ACDM growth history. 

The bulk flow enhancement ([I]) grows with increasing 
.D, since more extra dimensions imply more scalar fields 
on the brane contributing to the gravitational attraction. 
It also grows with decreasing r c , since smaller r c implies 
that departures from standard AD gravity turn on at ear- 
lier times. Because bulk flows measure today's evolution, 
higher- D models generate bigger departures from ACDM 
for fixed late-time normalization, since the force enhance- 
ment grows with D. 
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FIG. 1: The 1-dimensional rms velocity as a function of scale 
for standard gravity (black line) and brane-induced gravity, 
derived from linear analysis. The dashed (blue), dash-dotted 
(red), dotted (green), and top dashed (orange) lines show the 
maximum velocity allowed in our model with D = 5, 6, 7, and 
10, respectively, for late-time power spectra with ag ~ 0.88, 
consistent with observational constraints. The expansion his- 
tory and primordial normalization are identical in all cases. 

The enhanced gravitational attraction in cascading 
gravity also boosts the amplitude of density perturba- 
tions, resulting in a larger cr 8 at late times versus stan- 
dard gravity for fixed primordial normalization. The late 
time (7 8 is constrained by various observables, such as 
cluster abundance (76] [77] , the weak lensing power spec- 
trum [75] and galaxy clustering [75]. We discuss these 
constraints in Sec. |IV| The tightest constraint comes from 
X-ray galaxy cluster counts from the ROSAT All-Sky 
Survey [76], which for ACDM growth history imposes 
cr 8 < 0.88 at the 95% C.L. We therefore require that the 
late-time amplitude in our models approximately satis- 
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fies this bound. For a given D, this gives a lower bound 
on the allowed r c . On one hand, this is conservative: 
because our additional scalar force turns off in high- 
density environments, like galaxy clusters, the amplitude 
of the late-time linear power spectrum, when compared 
with data using ACDM methods, overestimates the clus- 
ter abundance [50]. On the other hand, the additional 
scalar force, if not completely screened, may increase 
the dynamical mass of any given cluster; this could in- 
crease the expected number of clusters at fixed dynam- 
ical mass, thereby tightening constraints. See [ST] for a 
careful study of dynamical effects in a variety of models. 
A similar treatment in our class of theories is work in 
progress [80] . 

Figure l] shows the 1-dimensional rms velocity vfL dim 
(the variance from zero mean of the theoretical distri- 
bution) as a function of scale for standard gravity and 
brane-induced/ Cascading Gravity with D = 5,6,7 and 
10, derived from linear theory using a Gaussian window 
function. For each D, we choose the minimum allowed 
value of r c , so each of the modified gravity curves have a 
similar late-time erg. 

To compare with data from the peculiar velocity sur- 
veys analyzed by Watkins et al. [T], we use a window 
function that includes observational effects but is only 
valid on 50/t^ 1 Mpc scales. Velocities obtained with this 
window function are denoted with superscript "W" . In 
practice, this results in a ~ 1% increase compared to the 
Gaussian window function, through the inclusion of some 
higher-fc modes. Over the range 5 < D < oo, we find 
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This is compared with the ACDM value, 179 km/s. The 
upper end of this range is nearly consistent at the 2a level 
with the observed bulk flow of 407 ± 81 km/s [J. While 
we have used linear theory to derive these results, we also 
perform N-body simulations of the non-linear evolution; 
these are in excellent agreement with the linear analysis 
for the scales of interest (c/. fjv]). 

The stronger effective gravitational attraction in our 
models also makes the occurrence of a high- velocity merg- 
ing system like the Bullet Cluster much more probable. 
Assuming that the majority of the infall velocity is caused 
by the ~ 1O 15 M main cluster, we treat the clusters as 
point particles released from rest from a large initial sep- 
aration (30 Mpc) and calculate the resulting velocity at 
5 Mpc, the initial separation for gas collision simulations. 
To include the rest of large scale structure, we add a bulk 
flow component to the infall velocity in quadrature. 

We find that the resulting velocity at 5 Mpc separation 
is 14% (D = 5) to 27% (D = 10) larger in our theories, 
where we again take the minimum allowed r c for each D. 
For D = 10 (and r c = 2750 Mpc), an initial velocity of 
3000 km/s at z — is a 4.8cr event, as compared with 
6.5cr in ACDM; or, in terms of probability, 6.6 x 10~ 7 
versus 3.3 x 10 -11 , a boost of over 10 4 . Meanwhile, at 
z = 0.5, closer to the actual merger redshift of z — 0.296, 
the modified gravity result is 3.9er from the mean, versus 



5.8ct in standard gravity. The probability of achieving 
that velocity is thus increased to 5.1 x 10~ 5 , compared 
to 3.6 x 10 -9 , again a ~ 10 4 probability enhancement. If 
the required initial velocity is reduced to 2000 km/s, the 
probability of such an occurrence in our model becomes ~ 
0.8% at z = and 14.2% at z = 0.5, respectively 257 and 
65 times more likely than the ACDM expectation (that 
is, respectively 2.4 and la in modified gravity, versus 4 
and 2.9cr in standard gravity). 



C. Modeling Cascading Cosmology 

Extracting exact cosmological predictions from a 
higher-dimensional set-up, such as Cascading Gravity, is 
technically very challenging; even the precise form of the 
modified Friedmann equation is not yet known. Hence 
our results must rely on some phenomenological input. 

Fortunately, the scales of interest for this study are 
well within the Newtonian regime. In this regime, the 
theory admits a local description on the 3-brane, which 
arises through a certain decoupling limit. In this limit, 
the complexities of the full higher dimensional theory 
that are irrelevant to AD physics are left out. The re- 
sult is an effective theory in AD with new degrees of 
freedom. The relevant degrees of freedom are weak-field 
gravity and D — A scalar fields coupled to the trace of 
the matter stress-tensor, describing brane-bending along 
each extra dimension. Whereas gravity is weakly cou- 
pled in this limit, the scalars have non-linear derivative 
interactions that are responsible for the Vainshtein ef- 
fect, a phenomenon generic to this kind of theory that 
leads to the approximate recovery of standard gravity in 
high density regions. While the precise form of the de- 
coupling theory for Cascading Gravity is not known, we 
draw upon known results in DGP to infer the non-linear 
interactions of the scalars. 

For the background, meanwhile, we assume a ACDM 
expansion history. As mentioned earlier, the form of 
the modified graviton propagator in Cascading Grav- 
ity suggests that brane-world corrections to the Fried- 
mann equation [83H86) are more slowly- varying functions 
of Hr c than in standard DGP [87]; they should closely 
resemble vacuum energy contributions. Furthermore, by 
assuming a ACDM expansion history, our analysis iso- 
lates the effects of the modified growth history [55] , We 
describe this approach in detail in Sec. [I] 

Finally, we wish to emphasize that, although our anal- 
ysis does not constitute a fully rigorous derivation of cos- 
mological predictions in Cascading Gravity, our results 
should most likely capture the essence of the predictions 
of the complete theory. Because the observables of inter- 
est are in the Newtonian regime, the description in terms 
of conformally-coupled scalar fields with non-linear inter- 
actions should be accurate. The details of the Vainshtein 
screening mechanism may vary, but the corrections to our 
predictions — which lie on scales where linear theory is 
valid — are expected to be small. In short, our approach 
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is not merely a toy model of cascading cosmology, but a 
first iteration in extracting predictions from a broad class 
of IR-modified gravity theories. 



I. MASSIVE/RESONANCE GRAVITY 

The defining feature of theories with infinite-volume 
extra dimensions, such as DGP and Cascading Gravity, 
is that AD gravity is mediated by a resonance graviton — 
a continuum of massive states — with general propagator 
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k 2 + m 2 (fc) 



(4) 



In DGP, for instance, m 2 (k) — r~ 1 k. In real space, this 
gives a gravitational potential that interpolates from the 
AD scaling, 1/r, at short distances to the 5D scaling, 
1/r 2 , at large distances, with a cross-over scale set by r c . 
In Cascading Gravity, the soft mass term m 2 (k) is a more 
complicated function of k, involving multiple cross-over 
scales [211 [22] . For simplicity, we shall assume that all 
cross-over scales are comparable and denote this common 
scale by r c . 

Because AD gravity is massive, each graviton has 5 po- 
larizations: the usual 2 helicity-2 states of GR, 2 helicity- 
1 states and 1 helicity-0 degree of freedom. At distances 
r <C r c , only the helicity-2 and helicity-0 degrees of free- 
dom are relevant — the helicity-1 states are very weakly 
coupled to matter and can be safely ignored. In the 
DGP model, the helicity-0 mode has a nice geometrical 
interpretation. It measures the extrinsic curvature of the 
brane in the extra dimension and is thus referred to as 
a brane-bending mode. Cascading Gravity theories have 
D — 5 additional hclicity-0 or scalar modes, accounting 
for the higher number of extra dimensions in which the 
brane can bend. This counting of degrees of freedom can 
alternatively be understood from a decomposition of the 
£>-dimensional massless spin-2 representation [2TJ [5S] . 

These scalar modes couple to the trace of the stress- 
energy tensor of matter on the brane and, combined with 
the hclicity-2 states, result in a one-graviton exchange 
amplitude between conserved sources having the tensor 
structure of D-dimensional massless gravity [2"2"] : 
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where we have neglected the graviton mass term for the 
scales of interest [55] ■ For non-relativistic sources, this 
corresponds to the modified Poisson equation 
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In other words, the gravitational attraction is a factor of 
1 + |j5f stronger than in standard Newtonian gravity. 

If© were valid in the solar system, the theory would 
already be ruled out by post-Newtonian constraints, for 



arbitrarily small m — this is the famous van Dam- 
Veltman-Zhakarov (VDVZ) discontinuity 90J . The res- 
olution, first proposed by Vainshtein for massive grav- 
ity ST] , is that the weak-ficld/linear assumption implicit 
in (5| is actually a poor approximation for the scalar 
modes in the vicinity of astrophysical bodies or in the 
early universe. Instead, as shown explicitly in the DGP 
model [921 193] . non-linearities in these modes in high- 
density regions result in them decoupling from matter, 
leading to an approximate recovery of Einstein gravity 
in, e.g., the solar system. 

Let us start with the DGP case. On scales -C H^ 1 
(and thus <C r c ), we can neglect time derivatives relative 
to gradients of the graviton helicity-0 mode. This mode, 
denoted by x, satisfies the approximate equation 



v 2 x - 
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where the V's denote spatial derivatives, and where 



/3 DGP = 1 + 2Hr c 1 



H 
3H 2 
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The overdot denotes a derivative with respect to proper 
time. The Vainshtein effect appears in two guises in the 
above. Firstly, in the early universe when Hr c 1, 
the coupling to matter density becomes vanishingly small 
since /?dgp 3> 1. Secondly, even at late times when 
Hr c < 1 and /?dgp ~ 1; sufficiently large overdensitics 
trigger non-linearities in \ and result in its decoupling. 

The analogue of Q is not yet known for Cascading 
Gravity. There are multiple scalars in this case, each of 
which is expected to exhibit Vainshtein screening. There 
have not yet been any successful calculations that keep 
the non-linearities of all scalar modes. For the purpose 
of this paper, we shall take a phenomenological approach 
and assume that all cascading scalar degrees of freedom 
obey an equation of the form ^ . This is consistent with 
our earlier assumption of a single cross-over scale r c . In 
other words, we collectively denote the scalars by x and 
assume that they satisfy 



V 2 x 



3/3a : 
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Using the scaling of the resulting deviation from GR in 
the solar system [59], we can infer that the cosmological 
deviation scales as j3 ~ H 2 r 2 . We then generalize ((8j) to 



ee 1 + 2H 2 r 2 c I 1 
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This modification makes a sharper transition in time 
than This is related to the fact that gravity rapidly 



5 



weakens on long length scales in models that can suc- 
cessfully degravitate the cosmological constant. This also 
means that whatever replaces ([9| in Cascading Gravity 
will likely exhibit even more efficient Vainshtcin screen- 
ing within collapsed objects. The effect of the x ne ld on 
matter is through its contribution to the actual potential 
"Jdyn that moves particles: 



*dyn = *N 
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where is the usual Newtonian potential. The In- 
dependence here is chosen to recover (JsJ) in the linearized 
limit. When performing N-body simulations, we solve ^ 
using a multi-grid relaxation scheme similar to the one 
described in [55] . Our N-body code is an updated version 
of the one used in previous work [30] , revised to solve ^ 
exactly and without resorting to the approximation of 
spherical symmetry. Meanwhile, we solve for in the 
usual way through Fourier transforms on a particle mesh. 
The code and further numerical results will be described 
in more detail elsewhere [50] , 

Most of our results are derived within linear theory, 
since this is a valid approximation for bulk flows over 
the scales of interest. In this regime, we can obtain a 
modified evolution equation for the density perturbations 
6 = Sp/p. Using the fact that the energy-momentum 
tensor on the brane is covariantly conserved, density and 
velocity perturbations evolve as usual via 



k k 
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These can be combined with the linearized version of Q 
to yield 
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Note that this is consistent with |6]) in the flat space limit 
/? — > 1. As another quick check, letting /? — > /?dgp an d 
setting D = 5, this agrees with explicit cosmological re- 
sults in DGP [93] . In this paper, we will solve for peculiar 
velocities using linear theory and cross-check the results 
with simulations. 



II. REVIEW OF BULK FLOW FORMALISM 

Since matter responds simply to gravitational gradi- 
ents at the linear level, it is straightforward to estimate 
the expected bulk flows on length scales where matter 
overdensities are in the linear regime. See §11.14 of [57] 
for an introductory exposition. 

The root-mean-square of the one-dimensional (1-dim) 
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on a scale R is given by 



P(k)W 2 (k 1 R)f 2 (k)dk 1 (14) 



where P(k) is the power spectrum of matter density fluc- 
tuations, and W(k) is a window function with scale R (To 
get the three-dimensional answer, we simply remove the 
1/3). Similarly to [55] , we choose the window function 
to be a Gaussian for the majority of our analysis: 



W^{k) = exp(-fc^) . 



(15) 



This is designed to capture only the small-fc/long- 
distance behavior of the power spectrum, i.e. the bulk 
flow. Note that the window function used for determining 
the bulk flow from peculiar velocity surveys [T] is slightly 
different from this, including some higher- kj smaller- R 
modes. We present this observational window function 
in Sec. |VI| and use it when comparing with the results 
of pp. In practice, however, it gives very similar results 
to the Gaussian window function. See Fig. [6] for a sneak 
preview of these window functions. 

The velocity also depends on the growth rate f{k) of 
perturbations, 



f(k) = dln #( a ' fc ) 
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where g is the growth function, and ao is the present scale 
factor. In the gravity theories of interest (as well as in 
theories with clustering dark energy) both g and / are 
generically scale-dependent. 

This formalism is similar to that used to set the nor- 
malization of the matter power spectrum. The parameter 
most commonly used for this purpose is a 8 , the matter 
fluctuation within 8ft. -1 Mpc spheres: 
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where W$(k) = 3j 1 (kR 8 )/kR 8 , with R 8 = 8ft -1 Mpc 
and ji a spherical Bessel function. Note that the window 
function for a 8 is a shallower function of k than that for 
a v . Hence, a 8 folds in more higher-fc modes, whereas 
peculiar velocities arc a direct probe of low-fc power. 

The long-distance modifications of gravity of interest 
boost peculiar velocities in two ways: 

• Faster development of structure at late times, en- 
coded in the growth rate /; 

• Greater build-up of structure due to the integrated 
influence of stronger gravity, resulting in a larger 
amplitude of the power spectrum (i.e., larger a 8 ). 

As mentioned earlier, however, the amplitude of the late- 
time matter power spectrum cannot be too drastically 
altered, as it is constrained by large-scale structure ob- 
servations. These will be taken into account in Sec. II VI to 
constrain our model parameters. Nevertheless, the boost 
in the bulk flow can be substantial because of the growth 
rate effect. In other words, whereas the power spectrum 
is the integrated result of the entire growth history, pe- 
culiar velocities are also sensitive to the present growth 
rate. 
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Bulk flows are often discussed in a variety of differ- 
ent ways, which can be confusing. The bulk flows that 
are measured using peculiar velocity surveys and other 
techniques are reconstructions of a full three-dimensional 
bulk flow of our region of space, but they are based on 
measurements that are separately only one-dimensional 
— our line-of-sight measurements always project the ve- 
locity of each object onto the radial direction. What 
is measured, then, is a 1-dimensional velocity for each 
object in the survey. By collecting and averaging over 
a large number of objects, however, we are able to re- 
construct the full 3D bulk flow [T], albeit with a small 
residual contribution from some smaller scales. This 
can be seen in the difference between our window func- 



tion ( 15 ) and the ideal experimental window function for 
surveys ( 32 1 , shown in Fig. [6j which includes a "bump" 
at around k ~ 0.05ft. Mpc -1 . 

Present observational techniques only allow reliable 
measurement of our local bulk flow. Thus, we have three 
independent measurements to compare with the expected 
variance in the local velocity calculated in ( fl4| ) . After re- 
viewing the results of our determination of the variance 
in Secs. |III|fV} we compare those theoretical expectations 
with the current best local bulk flow data in Sec. I VII Al- 
though it is an effect on a different length scale, we also 
discuss how our model change expectations for Bullet 



cluster initial velocities in light of [12] in Sec. VI 



III. ANALYTICAL ESTIMATES 

In this Section, we estimate the expected bulk flows 
in Cascading Gravity, first through analytical methods 
assuming matter-dominated cosmology (Sec. Ill A) and 



through fitting formulas of linear theory integration for 
ACDM expansion history (Sec.|IIIB). 



A. Enhancement in Matter-Dominated Cosmology 

In the linear regime and on scales much smaller than 
r c , density perturbations evolve according to ( 13 1. To get 



a rough estimate of the expected bulk flow enhancement 
analytically, we first ignore dark energy and consider an 
Einstein-de Sitter (f2 m = 1) universe. 

At early times (Hr c 1), the cosmological Vain- 
shtein effect results in f3 3> 1, and density perturbations 
evolve as in standard gravity, with growing mode solu- 
tion 8 ~ a. Once Hr c < 1, however, the extra scalar 
modes of Cascading Gravity become effective and (3 sa 1. 
Equation ( 13 ) then reduces to 
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with growing mode solution 



(19) 



where 7 is given by |2]). As a quick check, since 7 = 2/3 
for D — 4, this reproduces the usual matter-dominated 
growth 5 ~ a. 

Since H ~ t^ 1 in matter-dominated cosmology, the ex- 
cess growth from the onset of the modified gravity phase 
(when H = l/r c ) until the present time is 



6 = S stdsrav x (# r c r 7+ s . 



(20) 



v = v std grav x 



-7+j 



Through the continuity equation, ( 19 ) translates into v ~ 
5 = "/8/t for the peculiar velocity. The enhancement 
relative to standard gravity is thus 

37 
2 

This expression neatly captures the two contributions 
to the bulk flow excess described in Sec. [TTJ the larger 
growth rate, through the 37/2 factor, and the boost in 
the power spectrum amplitude through ( 20 1 . The above 



(H r c ) 



(21) 



derivation, which crudely assumes a sharp onset of the 
modified growth when H = r~ , actually agrees to within 
a few % with the exact integration of ( 13 ) for matter- 



dominated cosmology. This illustrates the rapid turn-on 
of the scalar force for Cascading Gravity, (10 1. 



B. Corrections from including dark energy 

The above analytic estimate, while useful in guiding 
our thinking, does not capture the weakening of growth 
triggered by the onset of cosmic acceleration. Using a 
ACDM expansion history in ( 13 ) is straightforward but 
requires numerical integration. The resulting growth 
rate, /(a, k), and power spectrum, -P(fc), are then sub- 
stituted back into ( 14 ) to obtain the velocity dispersion 



ct i] as a function of scale R. In comparing modified and 
standard gravity, we keep other cosmological parameters 
fixed at their fiducial values: il m = 0.24, il^ = 0.042, 
h = 0.73, n s = 0.96, and a primordial amplitude chosen 
to yield erg = 0.80 for ACDM growth history, consistent 
with WMAP 7-year data. 

In a similar spirit to Peebles' famous expression for 
the growth rate in the presence of dark energy, /acdm — 
OjJJ ■ , we use the results of numerically solving the linear 
perturbation equations to derive a fitting formula for the 
effect of our modifications to gravity on the amplitude of 
clustering today. For D = 5, ... 10 and over the range 
500 < r c < 7000 Mpc, we find that a reasonable fit is 
given by 



^8 



std grav , , 
Co x 



2 \Q-n(-7+!) 

- 3 H r c 



(22) 



We can also use a fitting technique to extend the expres- 
sion for the peculiar velocity on 50/i -1 Mpc scales given 
by to a formula valid over a range of scales ]J«r c : 



-a-dim 



(R) ~ 102 x f (§r c tf ) 



' 50h.~^Mpc 
, R [h-lMpc] 



l°-x 



(fg^l) 4 km/s,(23) 



7 



with 7 defined in 



IV. CONSTRAINTS ON THE AMPLITUDE 



The boost in the amplitude of density perturbations 
translates into a larger value for cr 8 compared to standard 
gravity, for fixed primordial amplitude. It is important to 
emphasize that the precise, early-time WMAP limit on 
a a applies only to the primordial amplitude as evolved 
to today using standard growth history; our choice of 
this amplitude, er 8 = 0.8, is consistent with WMAP 7- 
year results [75]- The amplitude of the late-time power 
spectrum is measured by various observations discussed 
below, constraining our models. 



D r c (in Mpc) 


5 


1665 


6 


2250 


7 


2486 


8 


2614 


9 


2694 


10 


2750 



TABLE I: Minimum values for r c in Mpc allowed by obser- 
vational constraints on the power spectrum amplitude. This 
contains the same information as Fig. [2] 



A. X-ray Clusters 

The tightest constraint comes from cluster counts using 
X-ray observations of the ROSAT All-Sky Survey [75] : 



cr,S 



Allowed region assuming amplitude enhancement < 10% 



0\24 



0.829 ± 0.0275 (Clusters) , (24) 




Number of dimensions 



FIG. 2: The shaded region of the D — r c parameter space 
shows the values allowed by observational constraints on the 
power spectrum amplitude. See also Table [I] 



Notable omissions in our discussion are large-scale 
tests, such as the cosmic microwave background temper- 
ature anisotropy and the Integrated Sachs- Wolfe galaxy 
cross-correlation. In DGP, these observables are the most 
constraining for r c [99 . However, analysis of these ef- 
fects requires evolving perturbations on horizon scales; 
this is beyond our Newtonian treatment. Though our 
decoupling-limit arguments (as in Sec. [i]) should give ro- 
bust predictions on Newtonian scales, horizon-scale ob- 
servables require a relativistic treatment of cosmological 
perturbation theory in Cascading Gravity. 



where the error bar includes a 9% systematic uncertainty 
in the mass calibration [75]. For our fiducial fl m = 0.24, 
this implies ag 0.88 at the 95% C.L. The translation 
from cluster abundance observations to as assumes stan- 
dard gravity, while in our theories the dependence of the 
mass function on the linear power spectrum amplitude is 
modified. We will come back to this point shortly. In this 
work we take the as constraint at face value and, with our 
fiducial choice <7g td grav = 0.8, impose as/af" d glav < 1.1. 
Using (22), this translates into a constraint on r c and D: 



o \°- 71 (-7+l) 



(25) 



The allowed region of the D—r c parameter space is shown 
in Fig. [2] Note that for a given number of space-time 
dimensions D, this translates into a lower bound on r c . 
In particular, we find r c > 1665 and 2750 Mpc for D = 5 
and 10, respectively. 

The above constraint is, on the one hand, conservative. 
Because our scalar force turns off in regions of high den- 
sity, such as galaxy clusters, we expect fewer and smaller 
non-linear structures in our model as compared with a 
ACDM model with identical present-day power spectrum 
normalization. On the other hand, the additional scalar 
force in our model leads to systematically higher dynami- 
cal masses in clusters |81j . This implies an increase in the 
number of clusters of fixed dynamical mass as compared 
with ACDM with the same late-time as- The bound 



we use ( 25 ) assumes that these opposing effects approxi- 
mately cancel. These considerations have been studied in 
detail by [55] for chameleon// (i?) cosmology. A similar 
treatment in Cascading Gravity is work in progress [80) . 

Substituting the minimum allowed value of r c for each 
D in ( 23 ) , we obtain the maximum allowed bulk flow as 
a function of scale. The result is shown in Fig. [I] for 



D = 5, 6, 7 and 10. On scales 50/i _1 Mpc probed by [J, 
using the appropriate observational window function (see 
Sec. VI I, we find the range Q: 



220 < v£_ db 



< 237 km/s . 



(26) 



This should be contrasted with the ACDM prediction of 
179 km/s. The upper range of (26) is therefore almost 



consistent at the 2a level with the observed bulk flow of 



407 ± 81 km/s pQ. In Sec. VI we present a more careful 
comparison with the data. 



B. Weak Lensing 

The Canadian France Hawaii Telescope weak Lensing 
Survey (CFHTLS) [78] finds 



erg 



(124 



0.5.S 



= 0.855 ± 0.086 (Weak Lensing) . (27) 



This is a much weaker constraint than (24 1 and is au- 



tomatically satisfied for the range of late-time values for 
as considered here. As a side remark, note that the con- 
straint on our modified gravity theories from weak lens- 
ing is actually weaker than implied by (27) because of 
the cosmological Vainshtein effect. The lensing kernel 
for CFHT sources spans the redshift range 0.25 < z < 1, 
where the effect of modifications to gravity is somewhat 
muted: for D = 10 and r c = 2750 Mpc, the case that 
yields the largest deviation from standard gravity, the 
difference from standard gravity varies from < 5% at 
z = 1 to «7.5% at z ~ 0.5. Hence the average enhance- 
ment of a% for observations made in this redshift range 
is less than the z = result suggests. 



C. Galaxy Clustering 

The amplitude of the late time matter power spectrum 
is also constrained by the clustering of galaxies measured 
by redshift surveys, such as the Sloan Digital Sky Survey 
(SDSS) PI!] and the 2 Degree Field Galaxy Redshift 
Survey (2dFGRS) [10 lj . However, an immediate com- 
parison of these observations is not straightforward due 
to the uncertainty in the galaxy bias. Over the mass 
range of interest, though, preliminary results using N- 
body simulations show little difference in halo bias in 
our model compared with the standard gravity predic- 
tion [50]. Marginalizing over bias, Seljak et al. [TpJ found 



cr 8 = 0.88 ± 0.06 (Galaxies) 



(28) 



Again this is a much less restrictive result than the X-ray 
cluster observations. It is worth noting that the value of 
a 8 that we use in our modified gravity results, which is 
in considerably better agreement with the peculiar ve- 
locity data than ACDM, coincides with the central value 
of ([28]). 



D. E g Parameter 

The expectation value of the ratio of galaxy-galaxy 
lensing to galaxy- velocity cross-correlations has been pro- 
posed as an observational test of gravity [102J. In linear 
theory, this combination, denoted by E g , is independent 
of bias and initial power spectrum normalization. Using 
SDSS luminous red galaxies at z = 0.32, [103] recently 
obtained 



Eg = 0.392 ± 0.065 . 



(29) 



In both ACDM cosmology and Cascading Gravity, this 
parameter is given by 



E 9 = 



fin 



/ ' 



(30) 



where fi m is the present-day matter density, and / — 
dlng/dlna is the growth rate at the redshift of observa- 
tion. With our fiducial value of fi m = 0.24, the ACDM 
result is E g — 0.387. Our IR-modified gravity theories 
predict 0.333 < E g < 0.346 for 10 > D > 5, assuming 
as before the minimum allowed r c for each D. Values in 



this range are consistent with ( 29 1 



Dependence on D, with fixed r c = 1000 Mpc 



o 



Q_ 

V> 

O 

o 

> 

■o 

£ 
o 

d) 
Q. 

X 




R in Mpc/h 



FIG. 3: The expected 1-dimensional peculiar velocity as a 
function of scale, using the Gaussian window function ( |15[ ). 
The curves are from linear theory whereas points (with 
bootstrap-estimated error bars) are from N-body simulations. 
The black curve/points are the standard gravity results. The 
long dashed (blue), dash-dotted (red) and dotted (green) 
curves/points are respectively the D = 5, 6 and 7 results, 
keeping r c — 1000 Mpc fixed. The points at a given scale 
are slightly offset from one another for readability For this 
comparison, <r 8 [ACDM] = 0.85. 



9 



Dependence on r c between 500 and 3000 Mpc, for fixed D=5 



— — — — — 3000 Mpc 

1500 Mpc 

■ 1000 Mpc 
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FIG. 4: The expected 1-dimensional peculiar velocity as a 
function of scale, using the Gaussian window function (15 1, 
for D = 5 and, from top to bottom, r c = 500, 1000, 1500 
and 3000 Mpc respectively in long-dashed (blue) , dash-dotted 
(red), dotted (green), and short dashed (orange) curves (linear 
theory), and circle, diamond, upright triangle, and inverted 
triangle points (simulations, with bootstrap-estimated error 
bars). The standard gravity results are plotted as a solid 
(black) curve and square points. The points at a given scale 
are slightly offset from one another for readability. For this 
comparison, <t 8 [ACDM] = 0.85. 



V. NUMERICAL SIMULATIONS 

While the linear theory analysis should be applicable 
to calculate bulk flows on scales > 50/i _1 Mpc, in this 
Section we explicitly check this against N-body simula- 
tions to include non-linear effects. Note that for the sole 
purpose of this comparison, we use a slightly higher pri- 
mordial normalization for which crg[ACDM] = 0.85. 

The N-body results presented here are obtained by in- 
tegrating the x equation of motion Q exactly, using a 
particle mesh approach. This is a notable improvement 
over earlier simulations presented in |30) . where a spher- 
ical approximation was employed. Although the exact 
evolution agrees with the approximate results reported 
in |30| to within < 5% for the power spectrum over 
the relevant scales, we nonetheless use the full code in 
the present work; a detailed comparison of the two ap- 
proaches will be discussed elsewhere [50] . 

We performed a series of particle-mesh simulations 
of 400/i _1 Mpc boxes on a 512 3 grid with 512 3 parti- 
cles. To determine bulk flow statistics, we perform a 
real-space average of velocities over a large number of 
spheres centered on points randomly placed throughout 



our simulation output, randomizing and rerunning each 
400/i _1 Mpc box 2 or 3 times for each choice of param- 
eters to gather better statistics. However, it is worth 
noting that each box gives highly consistent results with 
each of the other boxes and with the ensemble. We also 
use identical initial conditions across different parame- 
ter choices to isolate the effects of the new gravitational 
physics from random fluctuations. 

For each sphere, we calculate a Gaussian-weighted av- 
erage velocity, 



(v(R)[one sphere]) = — exp [-(r. t /i?) 2 ] 



(31) 



where is the distance from the ith point from the 
randomly-selected center, and N = ex P [ — ( r i/R) 2 ] 
is a normalization factor. This gives 3 independent 1- 
dimensional velocities for each sphere. All together, we 
average over 1200 — 4500 different velocities for each pa- 
rameter pair. We then calculate the standard deviation 
of all the sphere-averaged 1-dimensional velocities to ob- 
tain the rms velocity. 

Figures [3] and [4] compare linear theory integration 
(lines) with N-body results (points). Figure [3] explores 
the sensitivity to the number of extra dimensions, keep- 
ing r c fixed at 1000 Mpc. (For this value of r c , only the 
D = 5 case satisfies the observational constraints dis- 
cussed in Sec. |IV| these plots are made only to illustrate 
the sensitivity to D.) Figure. [4] shows the dependence on 
r c , fixing the number of dimensions at D — 5. In each 
case, we performed the calculations for R = 10, 25, and 
50/i _1 Mpc. We generated the error bars by bootstrap 
resampling subsets of the measured velocities, computing 
the rms for each subset, and using this as a dataset for 
constructing an error estimate. 

The upshot is that these figures show excellent agree- 
ment between linear theory and N-body simulations over 
the scales of interest. For bulk flows, the effects of non- 
linearities are basically absent, as hoped. This is in 
contrast with other velocity-related phenomena studied 
in |104j . where the effects of non-linearities persist to 
scales ~ 50/i _1 Mpc. 



VI. COMPARISON WITH DATA 

A. Bulk Flows 

Figure [5] compares the range of expected 1-dimensional 
velocity variances with the three local flow components 
reported in [V, . Our task in this Section is to quantify the 
extent to which these data are more likely in Cascading 
Gravity as compared to GR. We focus on the largest 
scales (50/i _1 Mpc) to minimize the influence of non- 
linear structures. Moreover, the data on different scales 
are correlated and would require a careful treatment of 
covariances. Note that all cosmological parameters are 
kept fixed in this analysis. We leave a comprehensive 
parameter likelihood analysis to future study. 
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As mentioned earlier, for the comparison with data on 
50ft, Mpc scales we use the same window function as 



determined in [1\ to analyze peculiar velocity surveys. 
This window function, shown in Fig. [6j is well-fitted by 



W 2 < w (k)~3- 



1.86 cos(93.5fc) + 0.0004 cosh(114.1 fc) - 11.58 sin(2fc) + 1.35sinh(46fc) 



4.89 cos(0.0012 k) + 0.73 cosh(121.3 k) 
I 



1.45 sinh(66.6fc) 



(32) 



Figure [6] compares this fitting function with the actual 
window function of 1 and with our Gaussian window 
function. (In Fi g.pl we instead use the Gaussian window 
function, since (|32|) only applies on 50ft -1 Mpc scales.) 



ties, which has vanishing mean and variance v\ 



vv 



-dim ' 



Allowed flow variance versus local data 
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FIG. 5: Comparison of the data from p] (points) with the 
range of bulk flows achievable in Cascading Gravity models 
(shaded region, using Gaussian window function 1 15}). The 
heavy curve defining the bottom of the shaded region is the 
result for standard gravity for our fiducial cosmology. Note 
that the 20/i -1 and 50/i -1 Mpc data points are not indepen- 
dent. Using only the 50h~ 1 Mpc data points and comparing 
their fit to standard gravity to the fit to the maximum allowed 
velocity (which is reached for D = 10, r c — 2750 Mpc), we 



find a Ax 



-4.6. 



The local bulk flow has three velocity components Vi, 
i G (1,2,3), measured by PQ. Their measurements have 
observational uncertainties cr,, which we assume repre- 
sent independent Gaussian error bars. Although these 
components are in truth related by a covariance ma- 
trix, for simplicity we assume that the data points are 
independent. This assumption has only a small effect 
on the statistics [105) . We compare these measurements 
with the expected distribution of 1-dimensional veloci- 



As- 

suming a Gaussian probability distribution, we can com- 
pare the different likelihoods through a simple x 2 statis- 
tic: 



i=l 



1-di) 



T 2 - 



(33) 



Recall that for our fiducial cosmology the standard 
gravity peculiar velocity variance is vj^^^ [std grav] = 
104 km/s. We wish to compare this with the maxi- 
mum allowed variance in our cascading gravity models, 
i>^ dim [mod grav] = 137 km/s, achieved with D = 10 and 
r c = 2750 Mpc. Comparing the x 2 values, we obtain 



2 

eff 



-4.6 



(34) 



in favor of the Cascading Gravity model, or —1.54 per 
degree of freedom. 

Although appreciable, this improvement in the fit does 
not by itself warrant the inclusion of two new parame- 
ters (r c and D) in a strict sense. However such parame- 
ter counting is often misleading. The underlying goal of 
long-distance modification of gravity is to relate these pa- 
rameters to an existing one, A. In a self-consistent anal- 
ysis of cosmological predictions in Cascading Gravity, it 
is conceivable that these parameters are in fact related to 
one another. For the purpose of this work, we can state 
conservatively that current peculiar velocity data tanta- 
lizingly hint at gravitational physics beyond ACDM but 
do not require it. 

Meanwhile, on the very large, 300ft -1 Mpc scales rele- 
vant to the result of [5] , our modifications to gravity can- 
not account for the observed bulk flow. For our D = 10 
and r c — 2750 Mpc model, we find up_ dim [mod grav] = 
29 km/s, versus 24 km/s for ACDM. 



B. Bullet Cluster 

As mentioned in the Introduction, recent hydronamical 
simulations of the Bullet Cluster have shown that an ini- 
tial velocity of v ~ 3000 km/s is required when the clus- 
ter and subcluster are ~ 5 Mpc apart to best reproduce 
X-ray observations. Lee and Komatsu [12] have recently 
estimated that the probability of such an initial velocity 
in the standard ACDM framework is between 3.3 x 10 -11 
and 3.6 x 10 -9 — respectively 6.5 and 5.8a from the mean 
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Comparison of Gaussian and Watkins et al. Window Functions 
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FIG. 6: Comparison of the Watkins et al. window function 
from [T] (squares), the fitting function given in ( |32[ ) (solid 
line) and the Gaussian window function (dashed line) with 
R = 50fe _1 Mpc, normalized such that their value is unity 
for k — > 0. The difference between the results for the two 
different window functions is very small, of order ~ 1%. 



of cluster velocities — , where the uncertainty comes from 
the evolution of velocities with redshift. 

Most of the infall velocity is caused by the gravita- 
tional attraction of the main cluster, which is estimated 
to have a mass of 10 15 M Q (That this is a good approx- 
imation for such massive clusters is established in |12|). 
We therefore compute the infall velocity in our models 
by treating the clusters as point particles released from 
rest from some large initial separation ~ 30 Mpc and in- 
tegrate their dynamics down to a separation of 5 Mpc as 
in 12J. For simplicity, we assume head-on collision. The 
gravitational potential due to the large cluster is obtained 
by solving ^ for \ and the standard Poisson equation for 
\1/n for a point particle of mass 1O 15 M . The motion of 
the subcluster is then obtained from the full gravitational 
potential given by (11). We also include in quadrature 



the effect of the enhanced bulk flow on 5 Mpc scales; 
however, the final velocity is almost entirely determined 
by the infall calculation. 

l ii 

We should expect the enhancement to scale as G^ s , 
where G c g is the effective Newton's constant enhanced 
by the additional scalar forces in our models — see ( 13 1. 



In the last stages of infall, however, the enhancement 
is somewhat tamed by the Vainshtcin mechanism. A re- 
lated point is that the enhanced force in modified gravity 
does not substantially reduce the velocity at 5 Mpc nec- 
essary to explain the merger velocity. Because of Vain- 



shtein screening, the final velocity of a particle falling 
from 5 Mpc to the center of a 1O 15 M cluster differs by 
at most 1% compared with the infall in standard gravity. 

We find that initial velocities are 14% to 27% larger 
in our framework, with the smallest difference coming 
from D = 5 and the largest for D = 10, again assum- 
ing the minimal value of r c allowed by constraints on a%: 
r c = 1665 and 2750 Mpc, respectively. These enhance- 
ments shift the mean of the probability density function 
for initial velocities calculated in [T2] , which is a function 
of logv. The larger bulk flow component at 5 Mpc also 
slightly widens the variance of the distribution. 

The end result is that the required initial velocities 
for the Bullet Cluster merger are much more likely in 
our model. Focusing on the D — 10, r c = 2750 Mpc 
case, the probability of an initial velocity of 3000 km/s 
at z = is increased to 6.6 x 10 -7 — 2. Ox 10 4 times more 
likely than the ACDM value; in other words, a change 
from a 6.5cr to a 4.8<r event. The probability of that 
initial velocity at z = 0.5, closer to the actual merger 
redshift of z = 0.296, is increased to 5.1 x 10~ 5 , again 
more than 10 4 times as probable as the ACDM result 
- a shift from 5.8c to 3.9o\ If we follow [12] and also 
consider the probability of finding an initial velocity of 
2000 km/s, the probability of such an occurrence in our 
model becomes 0.8% at z = and 14.2% at z = 0.5, 
respectively 257 and 65 times more likely than the ACDM 
expectation (that is, respectively 2.4 and la in modified 
gravity, versus 4 and 2.9<r in standard gravity). This 
significant increase in probability is another tantalizing 
hint of new gravitational physics. 



VII. CONCLUSIONS 

In this paper, we have explored how peculiar velocities 
are affected in a broad class of IR-modified gravity theo- 
ries called brane-induced gravity. On the scales of inter- 
est, these theories admit a local 4D description in terms 
of weak-field gravity plus D — 4 scalar fields coupled to 
the trace of the matter stress-tensor. These scalar de- 
grees of freedom effectively strengthen the gravitational 
attraction at late times and speed up structure forma- 
tion. As a result, peculiar velocities are systematically 
larger than those expected in standard gravity. Com- 
parisons between N-body simulations and linear theory 
calculations show that linear theory gives an excellent 
description of the physics of bulk flows. 

We have found that large-scale bulk flows can be en- 
hanced up to ~ 40% relative to ACDM cosmology. The 
enhancement is limited by observational constraints on 
us, the tightest limit coming from X-ray cluster abun- 
dance. The predicted peculiar velocities alleviate the 
tension with recent observations of the bulk flow on 
50/i _1 Mpc scales by Watkins et al. [U[2], from a ps 3<t 
discrepancy in ACDM gravity to a rj 2a difference in 
Cascading Gravity. The agreement between theory and 
data is improved by A\ 2 — —4.6 in our model. Although 
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modest, this improvement offers further motivation for 
more accurate bulk flow observations. Peculiar velocities 
are also enhanced on much larger scales (~ 300/i _1 Mpc) 
probed by kinetic Sunyaev-Zel'dovich observations, but 
this is insufficient to explain the enormous bulk flow in- 
ferred by OS]. 

Violent merging systems, such as the Bullet Cluster, 
are much more probable in Cascading Gravity. Drawing 
on the recent analysis of [12], we have found that the 
occurrence of the Bullet Cluster in our theories is > 10 4 
times more likely than in standard gravity. 

We are pursuing various parallel tracks to improve 
upon the preliminary analysis presented here. As men- 
tioned above, the tightest constraint on our models comes 
from cluster counts. Through N-body simulations and 
various semi-analytical techniques, we are currently de- 
termining the halo mass function in Cascading Gravity 
as a function of D and r c [80]. This will allow a more 
accurate comparison with X-ray data. With regards to 
the Bullet Cluster, we can use similar simulations to cal- 
culate the distribution of initial velocities for subclusters, 
following the standard gravity analysis of [T2"] . 

On the theoretical side, the calculations presented here 
rely on some phenomenological input. While we are con- 
fident that our results capture the essence of Cascading 
cosmology predictions, they do not derive from a rigorous 
cosmological analysis of the complete higher-dimensional 
theory. An important first step in this direction would 
be to obtain the decoupling limit of Cascading Gravity 
while keeping the non-linearities in all scalar degrees of 
freedom. This would allow a derivation of the modified 
Friedmann equation, as in [3T] for DGP, as well as the 
perturbation equations in the Newtonian regime. This is 



work in progress [106) . 

The observations of large bulk flows and violent clus- 
ter mergers offer tantalizing evidence that structure is 
evolving more rapidly than predicted by ACDM cosmol- 
ogy. As we extend our measurements of large-scale bulk 
flows beyond our local region and discover an increasing 
number of merging systems, it will become clear whether 
these are statistical flukes or the first indication of a new 
realm of gravitational physics on cosmological distances 
in the late universe. 
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